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Abstract 

• We examine quenched chiral logarithms in lattice QCD with overlap 

Dirac quark. For 100 gauge configurations generated with the Wilson 
K/i \ gauge action at /3 = 5.8 on the 8^ x 24 lattice, we compute quenched 

^ • quark propagators for 12 bare quark masses. The pion decay constant is 



extracted from the pion propagator, and from which the lattice spacing 
is determined to be 0.147 fm. The presence of quenched chiral loga- 
rithm in the pion mass is confirmed, and its coefficient is determined 
to be (5 = 0.203 it 0.014, in agreement with the theoretical estimate in 
quenched chiral perturbation theory. Further, we obtain the topological 
susceptibility of these 100 gauge configurations by measuring the index 
of the overlap Dirac operator. Using a formula due to exact chiral sym- 
metry, we obtain the r/' mass in quenched chiral perturbation theory, 
rrir,' = (901 ± 64) MeV, and an estimate of 5 = 0.197 ± 0.027, which is 
in good agreement with that determined from the pion mass. 
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1 Introduction 



In quenched chiral perturbation theory @], it has been asserted that some 
quantities in QCD possess logarithmic dependence on the bare quark masses 
as the latter approach zero, similar to the unquenched case. In particular, 
for the pseudoscalar meson mass and the condensate, there are extra chiral 
logarithms due to the r/' loop void of topological screening in the quenched 
approximation. Explicitly, the pion mass to one-loop order can be written as 

ml = Cmjl - 5[ln(Cm,/Aj) + I]} + Bm] , (1) 

where denotes the bare ( u and d ) quark mass, the chiral cutoff which 
can be taken to be 2A/27r/^ ( ~ 132 MeV ), C and B are parameters, and 
5 the coefficient of the quenched chiral logarithm. 
Theoretically, 5 can be estimated to be 



ml, 



(2) 



where m^' denotes the t]' mass in quenched chiral perturbation theory, /,r the 
pion decay constant, and Nf the number of light quark flavors. For m^/ = 900 
MeV, = 132 MeV, and Nj = 3, (D gives 6 ~ 0.2. Evidently, if one can 
extract 6 from the data of in lattice QCD, then m^/ in quenched chiral 
perturbation theory can be determined by (0). 

Presumably, the chiral logarithm in (|I]) can be observed unambiguously in 
quenched lattice QCD. However, it is rather difficult to disentangle a logarithm 
from a power series. So far, there is no compiling evidence to single out (|I]) 
among many functional forms which seemingly can fit the data of very well 
( e.g., ml = A + Cmg + Bm^ ). 

For lattice QCD with Wilson-Dirac quark, the chiral symmetry is explicitly 
broken and the quark mass is additively renormalized, thus one does not expect 
that quenched chiral logarithm can be identified unambiguously at finite lattice 
spacing. In fact, if one tries to fit (0) to the data of with Wilson-Dirac 
quark, one would obtain the coefficient of quenched chiral logarithm, 5 ~ 
0.06 — 0.12 0, quite smaller than its theoretical expectation 0.2. For lattice 
QCD with staggered quark, the coefficient of chiral logarithm also turns out to 
be rather small, 5 ~ 0.06 [^, similar to the case of Wilson-Dirac quark. These 
discrepancies indicate that quenched chiral logarithm may not be properly 
reproduced in lattice QCD with these two lattice fermion schemes, at finite 
lattice spacing. 

With the realization of exact chiral symmetry on the lattice, one expects 
that quenched chiral logarithm can be reproduced in lattice QCD with overlap 
Dirac quarks 0, 0|. In a recent study p|, the coefficient of quenched chiral 
logarithm 5 is estimated to be in the range 0.15 — 0.4 for the chiral cutoff 
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in the range 0.6 — 1.0 GeV, by fitting (|1|) to their data of m^. At K^a = 
2\/27r/^a = 1.17 GeV, their data seems to give 6 ~ 0.41(9), which is too large 
comparing with the theoretical estimate, 5 ~ 0.18 — 0.20. Thus, it is not clear 
whether quenched chiral logarithm has been unambiguously identified. 

Besides from the data of m^, one can also obtain S by extracting m^/ 
from the propagator of the disconnected hairpin diagram. However, to com- 
pute the propagator of hairpin is very tedious. Fortunately, due to the chiral 
symmetry of Dc (^) in the quark propagator {Dc + mq)~^, only the zero modes 
of Dc can contribute to the hairpin diagram, regardless of the bare quark mass 
rriq. Thus one can derive an exact relation (|l9|) between the r]' mass in quenched 
chiral perturbation theory and the index susceptibility of any Ginsparg- Wilson 
lattice Dirac operator, without computing the hairpin diagram explicitly. Ex- 
plicitly, it reads as 



< = (3) 

where Nf denotes the number of light quark flavors, V the space-time volume, 
and — n^Y) jV the index susceptibility of any Ginsparg- Wilson lattice 
Dirac operator in the quenched approximation. Then @ and @ together 
gives 



{{n+-n_Y) 



(4) 



where Ng is the total number of sites of the lattice. A salient feature of is 
that S can be determined at finite lattice spacing a, by measuring the index ( 
susceptibility ) of the overlap Dirac operator, and with f^^a extracted from the 
pion propagator. Further, we suspect that the value of 6 (^) is scale invariant 
at least for a range of lattice spacings including the continuum limit a — 0. 

Now it is clear that, in order to confirm the presence of quenched chiral 
logarithm in lattice QCD, one needs to check whether the coefficient 6 de- 
termined from the data of agrees with that (^) obtained from the index 
susceptibility. Further, they should be close to the theoretical expectation 
S ~ 0.2. This is a requirement for the consistency of the theory, since the 
quenched chiral logarithm in (|l|) is due to the rj' loop coupling to the pion 
propagator through the mass term ( in the chiral lagrangian ), thus S ( m^/ ) 
must be the same in both cases. Otherwise, a consistent picture of the theory 
has not been established, and one may infer that something must have gone 
wrong, either in the formulation or in the implementation of the theory. We 
regard this consistency requirement as a basic criterion for lattice QCD ( with 
any fermion scheme ) to realize QCD chiral dynamics in continuum. 

In this paper, we examine whether lattice QCD with overlap Dirac quark 
can satisfy this basic criterion. 
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2 Pion mass and decay constant 



In this section, we set up our notations for computing quenched quark propa- 
gator as well as meson propagators in lattice QCD with exact chiral symme- 
try. Then we describe the practical implementation for overlap Dirac quark 
in our computations. Finally, we present our results of m^a^ and /ttCi for 12 
bare quark masses, and determine the coefficient of quenched chiral logarithm, 
S = 0.203 ±0.014. 



2.1 Some general formulas 

First, we recall some general formulas in lattice QCD with exact chiral symme- 
try. For any massless lattice Dirac operator D satisfying the Ginsparg- Wilson 
relation [^] 

D^, + ^,D = 2raD^,D , (5) 
it can be written in terms of a chirally symmetric Dirac operator Dc ||10|| , 

D = D,{l + raD,)-' . (6) 
Then the bare quark mass is naturally added to the Dc in the numerator of 

Dimg) = {Dc + mg)(l + raDc)-' . (7) 
Thus the quark propagator is 

{Dc + m,)-i = (1 - rmga)-'[D-\mg) - ra] , (8) 

which is chirally symmetric in the massless limit. This exact chiral symmetry 
is the crucial feature one would like to preserve for any quark coupling to 



physical hadrons [11 



Then the pion propagator can be written as 

M{x,t- 0,0) = ti{j,{Dc + mg)-\0,x)-f,{Dc + mg)-\x,0)} , (9) 

where the trace tr runs over Dirac and color space. Note that any one of the 
quark propagators in (|^) can form the rj' loop which leads to the quenched 
chiral logarithm in the pion mass as well as the condensate. 

With the pion propagators, one can compute its time correlation function 

G{t) = Y.M{x,t;6,0) , (10) 

which can then be fitted by the usual formula 

G^{t) = ^_[e-™-'^* + e-"^-'^(^-*)] , (11) 
2m7ra 
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to extract pion mass m^a, and pion decay constant, 



f^a = 2mga-!-^ . (12) 



mta? 



Besides the exact chiral symmetry, we also require that iDc is hermitian, 
similar to the Dirac operator 27^ (9^ + iA^) in continuum. Then we have 
D\ = 75Z}c75- As a consequence of hermiticity and chiral symmetry, both 
quark propagators in (H) can be expressed in terms of that propagating from 
(0, 0) to (x, t), i.e., 

M(f,t;0,0) = ti{[iD, + mg)-\0,x)]\D, + m,)-\x,0)} 

= E j: [{D. + rn,)-Z{x,0)y{D^ + m,)-Z{x,0) (13) 

a, 6=1 a, 13=1 

where a, b denote color indices, and a, (3 Dirac indices. In other words, to 
obtain the time correlation function (0) of pion, one only needs to compute 
12 columns of the matrix [D^ + mq)~^. 

2.2 Practical implementation for overlap Dirac quark 

The massless overlap Dirac operator reads as 

D = moa-Ml + 75-^| (14) 



where denotes the hermitian Wilson-Dirac operator with a negative pa- 
rameter —mo, 

= 'j^D^ = 75 (-mo + 'j^t^ + W) , (15) 



7^t^ the naive fermion operator, and W the Wilson term. Then D (|14|) satisfies 
the Ginsparg- Wilson relation (^ with r = l/(2mo). Note that the value of mo 
has to be chosen properly such that the overlap Dirac operator can capture the 
topology of the gauge configuration, as well as behaving like a massless Dirac 
fermion even in a topologically trivial gauge background. In other words, on 
a finite lattice, the proper range of mo for any gauge configuration is smaller 
than that (0 < mo < 2) in the free fermion limit. However, one usually does 
not need to finely tune mo except for rough gauge configurations at strong 



couplings ||T2[. In this paper, we fix mo = 1.3 for our computations. 

Now the difficult task is to compute the quark propagator with exact chiral 
symmetry (§), which amounts to computing 
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where the positive hermitian operator D{'mq)D\mq) can be written as 



D{mq)D\mq) 



m„ 



<+{<--r 2 + 75- 



+ ^=75 



(16) 



Then the quark propagators from the reference point (0, 0) to all sites (x, t) 
can be obtained by solving 12 column vectors Yca{x^t) ( where c and a are 
color and Dirac indices of the quark field at (0, 0) ) in the following linear 
system for 12 point source vectors on the r.h.s., 

D{mq)D\m,)Y,^{x, t) = , (17) 

where D{mq)D'^{mq) is expressed by (|16|). Due to the rather small size of 
the physical memory in the present generation of computers, the viable ways 
to solve the linear system ([T7| ) are iterative methods, in which the conjugate 
gradient algorithm is the most optimal for positive definite matrices. Starting 
from an inital guess of Y^^\ we multiply DD"^ to Y^^\ and update y^^-* to Y^^^ 
according to the conjugate gradient algorithm, then iterate this process until 
converges to the solution with desired accuracy e, i.e., 

\\D{mq)D\mq)Y}:\x,t) - < e . (18) 

In this paper, we set e = 10^^^. 

Now the updating process involves the multiplication of the matrix 

1 



(19) 



to the vector Y^^\ However, (19) does not have analytic closed form. This 
is the major challenge for lattice QCD with overlap Dirac quark. A way to 
proceed is to express (|19D in terms of a rational approximation |l! 



E 



(20) 



where hi and di are some positive definite coefficients. Then the matrix-vector 
multiplication 



1 



1 = 1 "T 1 = 1 



can be evaluated by invoking another conjugate gradient process to the linear 
systems 



{Hl + di)Zi = Y^^, 1 = 1, 



k. 



(21) 
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Instead of solving each Z[ individually, one can use multi-shift CG algo- 



rithm [|I4|, and obtain all Z[ altogether, with only a small fraction of the total 
time what one had computed each Zi separately. Evidently, one can also ap- 
ply multi-shift CG algorithm to (0) to obtain several quark propagators with 
different bare quark masses. 

Now the computation of overlap Dirac quark propagator involves two nested 
conjugate gradient loops : the so-called inner CG loop (^TD , and the outer CG 
loop ([T7|). The inner CG loop is the price what one pays for preserving the 
exact chiral symmetry at finite lattice spacing. 

The computational cost can be saved more than 50% by observing that 
D{vrLq)D\mq) (|16D commutes with 75. Thus one can choose chiral sources on 
the r.h.s. of fp] 



X, t) 



where p{x, t) = ^^St^ for a point source. Then we have 



75 Ic 
75 Ic 



a 



-lea, a 



1,2 
3,4 



and 



75 ^« 



a 
a 



1,2 
3,4 



Now one of the two matrix- vector multiplications in the outer CG loop can be 
eliminated, since the l.h.s. of (ITTD becomes 



ml + 



2ml 



1 



Yea (22) 



where '+' for a = 1,2, and '-' for a = 3,4. Further, due to the projector 
(75 ± l)/2, half of the matrix- vector multiplication in 



H,, 



--Yr, 



can be saved. 

The next crucial question is how to fix the values of coefficients bi and 
di in ( pOD such that the inverse square root function can be approximated 
optimally. Fortunately, this problem has been solved by Zolotarev [|l^ in 1877, 
using Jacobian elliptic functions. A comparative study of Zolotarev's optimal 
rational approximation versus other schemes for the overlap Dirac operator 



has been reported in Ref. [16 . 
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For the inverse square root function 

1 



X 



X 



^ [1) ^max] ; (23) 



the Zolotarev optimal rational approximation is 



(24) 

Ui=i{x + C2I-1) 



where 

sn\lK/2k; k 



ci = 



— 1 — 1 /Xmax ) (25) 



l-sn'^{lK/2k;K) 

the Jacobian elliptic function sn(M, k) = r] is defined by the elliptic integral 

u{v) = r , (26) 



K = m(1) is the complete elliptic integral, and do is uniquely determined by 
the condition 

max[l - \/x/(x)]|^e[i,^_,] = -min[l - Vxf{x)]\^^[i,^^^^] 

For our purpose, it suffices to fix do such that /(I) = 1. 
Since the eigenvalues of Hy^ are bounded, 

Amm ^ |A(i/u,)| < Xmax , (27) 



we can rescale H^, 



hyj — Hyjj Xfnin 5 (28) 



such that the eigenvalues of /i^ fall inside the interval [1, {Xmax/ Xmin)^]- 

After obtaining the partial fraction of (|2^) , we have the optimal rational 
approximation to the inverse square root of 



1 A bi 



IH^^ Xmin 1=1 hi, + d 

where 



(29) 



di = C21-1 , (30) 

, , I\i=l {C2i - C2I-1) , s 

' = "^Off^ fr r ^ ' 

lli=l,i^l[C2i-l - C2I-1) 
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and the parameter k in the Jacobian eUiptic function sn[lK/2k; n) is 



K = \/l - (Amin/AmaaO^ ■ (32) 

Note that, for a given order k, the smaller the interval [1, {Xmax/ ^min)^], 
the more accurate the rational approximation is. Thus, it is advantageous 
to narrow the interval [Xmin, Xmax] by projecting out the eigenmodes at both 
ends. Since the spectrum of is very dense near the upper bound, we only 
need to project out one or two eigenmodes near Xmax, just to get the precise 
value of Xmax- On the other hand, we usually project out 20 or more low-lying 
eigenmodes of H^, such that {Xmax/ Xmin)^ is always less than 2500. Then we 
project the inner CG loop to the complement of the vector space spanned by 
these eigenmodes. This not only gives a better rational approximation for the 
inverse square root of H^, but also reduces the number of iterations of the 
inner CG loop, which is proportional to the condition number {Xmax / Xmin)^ ■ 

We use Arnoldi algorithm to compute a selected subset of eigenmodes of 
H^, which correspond to high and low lying ones of H^. Denoting these 
eigenmodes by 

H^Uj = XjUj, J = 1, ■ ■ ■ , , (33) 

we project the inner CG loop (^) to the complement of the vector space 
spanned by these eigenmodes ( i.e., multiplying both sides of (^) by the 
projector P = 1 — J2]=i UjV^p and use PH^ = H^P ), 

n 

{hi + di)Z^ = F« = (1 - E ■ (34) 

i=i 

Note that Xmin and Xmax in Eqs. (|27| ) - (|32|) now refer to upper and lower bounds 
of \X{Hu,)\, where 

n 

Hyj = PH^ = H^-Y^ XjUju] . (35) 
i=i 

Then the matrix-vector multiplication at each step of outer CG loop can be 
expressed in terms of the projected eigenmodes (|33|) plus the solution obtained 
from the inner CG loop ( p^ in the complementary vector space, i.e.. 



1 1 _ _ " A 

i/^^yw ^ -^H^ Yl biZi + E -nh^^^lY^"^ = S'^'^ ■ (36) 

The accuracy of S^^^ can be measured by the deviation 

y(i)ty(i) 
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which would be zero if ( |36D is exact. 

In this paper, we project out 20 lowest-lying and 4 largest eigenmodes of 
H^, and use 20 terms ( = 20 ) in the Zolotarev optimal rational approxima- 
tion, and set the convergence criterion for inner CG loop to 10~^^, then cTj is 
always less than 10~^^ for any iteration of the outer CG loop. 



2.3 Results 

With the Wilson SU{3) gauge action at /? = 5.8, we generate 100 gauge con- 
figurations on the 8^ x 24 lattice, using Creutz-Cabibbo-Marinari heat bath 
algorithm p3|, For each configuration, we compute quenched quark prop- 



agators for 12 bare quark masses. Then the pion propagator (|T^ and its time 
correlation functions ([T0|) are obtained, and the latter are fitted by ([TT|) to 
yield the pion mass mT^a and decay constant /jr^- 

In Figs. ID and|^, we plot the pion mass square [rriT^a)'^ and decay constant 
/^a versus the bare quark mass niga respectively. The data of /^a can be fitted 

by 

f^a = 0.0984(3) + 0.1635(15)(mga) (38) 

with x^/d-O-f- = 0.03. Now we take the value f-^a = 0.0984(3) at m^a = 
equal to the experimental value of pion decay constant f-j^ = 132 MeV times 
the lattice spacing a, 

Ua = 0.0984(3) = 132 MeV x a , (39) 

then this gives an estimate of the lattice spacing 

a = 0.147(1) fm . (40) 

Thus the lattice size is ~ (1.2 fm )^ x 3.5 fm . Since the smallest pion mass 
is ~ 418 MeV, the lattice size is ~ (2.5)^ x 7.5, in units of the Compton 
wavelength ( ~ 0.47 fm ) of the lowest-mass pion. 

Fixing the cutoff = 2-\/2vr/7r as in quenched chiral perturbation theory 
( i.e., A^a = 2v^7r/^a = 0.8745 ), we fit (|l]) to our data of (m^a)^, and obtain 

S = 0.2034(140) , (41) 
Ca = 1.1932(182) , (42) 
B = 1.1518(556) , (43) 

with x^/d.o.f. = 0.03. 

Even though one may not easily detect the quenched chiral logarithm in 
Fig. [D, it can be unveiled by plotting {mT.aY' / {niqa) versus m^a, as shown in 
Fig. ^ Further, we can subtract the quadratic term Brn?^ [ with the value 
of B given in (^3|) ], and check its dependence on In(mg) explicitly. In Fig 
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we plot {m.„aY / {niqa) — BijUqa) versus log(mqa). The presence of quenched 
chiral logarithm is evident. 

In passing, we note that a finite lattice must impose a lower bound for 
the pion mass, as well as the corresponding one for the bare quark mass. If 
one decreases the bare quark mass beyond its lower bound, then the resulting 
pion mass would be fiattening ( or increasing ) rather than decreasing. This is 
essentially due to the finite volume effects of the zero modes of the overlap Dirac 
operator in the pion propagator, which is proportional to |n+ — n^\/ {m^a^Ns) 
as TUqa — 0. Thus, only in the infinite volume limit ( A^^ ^ oo ), one can 
obtain zero pion mass with zero bare quark mass. Nevertheless, for a finite 
lattice, the finite volume effects of the zero modes can be suppressed if the bare 
quark mass m^a ^ \J\''^+ ~ "^-l/^s- For a lattice of size 8^ x 24 at /? = 5.8, 
we find that m^a > 0.06 is sufficient to suppress the finite volume effects of 
the zero modes. This can be verified by comparing the results ( 
those obtained in larger volumes, e.g., 12^ x 24 at = 5.8 ( which is about 
3.4 times of the present volume ). Another plausible check is to see whether 
the 5 extracted from the pion mass agrees with that obtained from the index 
susceptibility, Eq.(4), since the index susceptibility and the /^a ( see Fig. 2 
) both presumably suffer the least from the finite volume effects. The latter 
check is performed in Section 3, and the agreement between these two 5's 
suggests that the finite volume effects may be under control. 

3 Index susceptibility and the rj' mass 

In this section, we derive a relationship between the index susceptibility of 
any Ginsparg- Wilson lattice Dirac operator and the r]' mass as depicted in 
quenched chiral perturbation theory. Then we measure the index of overlap 
Dirac operator for those 100 gauge configurations used for computing quark 
propagators in Section 2, and obtain the index susceptibility, the t]' mass, 
and an estimate of 5 = 0.197 ± 0.027, which is in good agreement with the 5 
determined from the pion mass in Section 2. 

3.1 A formula due to exact chiral symmetry 

The propagator of ( fiavor-singlet ) r]' can be written as 

-tr{75(D, + mg)-i(|/,y)}tr{75(D, + m,)-i(x,x)} (44) 

where the second term on r.h.s. is the commonly called hairpin diagram. Note 
that the first term is just the usual pion propagator (^, while the second 
term can be regarded as pion propagator with the t]' mass insertion due to the 



41)- 



to 



10 



nontrivial interactions of gluons between two seemingly disconnected quark 
propagators. 

In the quenched approximation, the expectation value of the hairpin dia- 
gram in momentum space can be written as 

E e*''^'"'^ tr{75(I^c + m,r\y, y)}ii{i,{D, + m,)-\x, x)}\ 
1 ^1 



-m^— -VZ (45) 

pj. _|_ / _|_ 

where the brackets (■ ■ ■) denote quenched average over gauge configurations 
with weight e""^" ( Ag : pure gauge action ), and 

2mq 



as defined in (pUf ) and (|12D . Equation (^Sf ) can be regarded as the definition of 
the rj' mass in quenched chiral perturbation theory. 
At zero momentum p = 0, (|45|) becomes 

(^E tr{75pc + m,)-i(2/,y)}tr{75pe + m,)-i(x,x)}^ = ^ (46) 

Now the expression inside the brackets (■ ■ ■) can be evaluated exactly in 
terms of the index of { D ). Explicitly, 

^ tr{75(Z^c + m,)-\x, x)} = Tt{j5{D, + m,)"^} 

X 

= (1 - rmgtt)"^ ^0)^750a{[mg + (1 - rmga)Aa]"^ - ra} (47) 

a 

where (|^) and (H) have been used, and the trace Tr is evaluated in terms of the 
eigenvalues {Xa} and eigenvectors {4>a} of the Ginsparg- Wilson lattice Dirac 
operator D (^) with Dc satisfying {iDcY = iDc- Recalling the properties of 
the eigensystem of D [p!7| . 



for Aq, 7^ A* 



±1 for Xa = A* 



one immediately sees that only the real eigenmodes A = and {ra)~^ can 
contribute to (^Tf). Since the factor {[m, + (1 — rmga)Xa]~^ — ra} vanishes for 
Xa = {ra)~^, thus only the zero modes contribute to (0), i.e.. 



Y,ti{^,iD, + m,)-\x,x)} = . (4. 
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Substituting ( ^HD into l.h.s. of (^61), one obtains 

iv/<"^ - "->'> = ^ , (49) 

where the factor Nf denotes the number of hght quark flavors, which accounts 
for the number of hairpin diagrams contributing to the r]' mass. Since ( ^Of ) 
is an exact result following from the definition of rrir]' in (145|) , it offers the 
best way to obtain the f]' mass in quenched chiral perturbation theory. One 
just measures the index ( susceptibility ) of the Ginsparg- Wilson lattice Dirac 
operator, without computing the hairpin diagram at all. This is one of the 
advantages of preserving exact chiral symmetry on the lattice. Note that (|49D 
is independent of the bare quark mass rriq, thus the 1]' can remain massive even 
in the limit rriq —>■ 0. This is the crucial feature of the rj' to distinguish itself 
from the octet of approximate Goldstone bosons, the vr, K and t]. 

Recently, (|49|) is derived pO[ differently starting from the anomalous fiavor- 



singlet axial Ward identity, which is formally equivalent to replacing the con- 



tinuum topological charge density p in the Witten-Veneziano formula |]18|, |l£ 



4Nf f 1 
"^V = 7/ y d''x{p{x)p{0)) , p = ^e^^uXatiiF^uFxa) , (50) 

with the axial anomaly q = tr[75(l — raD)] of the overlap Dirac operator, 
and the integration / d'^x with the summation over all sites J2x- That is, the 
integral in (|50D is formally transcribed as 

/ d'x{pix)pm ^ a-' E(5(^)5(0)) ■ (51) 

Since (^0[) is derived in the large Nc limit ( which is valid to lowest order in 
1/Nc ), a relevant question is what value of Nc in Eq. (|50|) is supposed to be 
valid. Further, even though the r.h.s. of (^Tj) is well-defined on a finite lattice. 



some subtleties ||2T| could emerge as the 2-point function {q{x)q{0)) ( x 7^ 
) becomes non-positive in the continuum limit, since p{x) is odd under time 
reflections. Thus, a real m^/ would require a divergent contact term (p^(0)). 
How to fix this contact term is the basic problem of Witten-Veneziano formula 
( pOD in the continuum. On the other hand, for a finite lattice, one can use the 
well-known relation J2x li^) = Uj^ — n- to turn the r.h.s. of (|5l|) into the index 
susceptibility, 

a-'j:{q{x)qm = ^EEi^i^Hv)) = ^(K -^-)') , (52) 

X ^ X y " 

which is well-defined in the limit V 00 and a 0. Nevertheless, it is still 
interesting to see whether (g(x)g(O)) ( x 7^ ) would become negative in the 
continuum limit, together with a divergent contact term. 
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It should be emphasized that our derivation of (^) is quite different from 
that of Ref . pOf . We start from Eq. (^51) which serves as the definition of the 
7]' mass in quenched chiral perturbation theory. Then the subsequent steps 
leading to (^Q]) are exact due to chiral symmetry, without any approximations 
or assumptions. This is in contrast to the approach used in Ref. 



m 



which the anomalous fiavor-singlet axial Ward identity is used, in the limit 
u = Nf/Nc — > 0. The subtleties arising from {q{x)q{0)) in the continuum limit 
cannot be resolved unless one transcribes the expression in terms of the index 
susceptibility of a Ginsparg- Wilson lattice Dirac operator. Finally we note 



that a derivation of (^9]) similar to our approach is also presented in Ref. ||22 
for the overlap Dirac operator. 



3.2 Measurement of the index of the overlap 



The index of overlap Dirac operator ([T^) is 



mdex{D) = ^tr[75(l - arD{x,x))] = hh^ - h+) (53) 

where h^{h^) is the number of positive ( negative ) eigenvalues of the hermi- 
tian Wilson-Dirac operator (p^. However, one does not need to obtain all 
eigenvalues of in order to know how many of them are positive or nega- 
tive. The idea is simple. Since has equal number of positive and negative 
eigenvalues for mo < 0, then one can just focus at those low-lying ( near zero ) 
eigenmodes of H^, and see whether any of them crosses zero from positive to 
negative, or vice versa, when mo is scanned from up to the value ( e.g., 1.30 
in this paper ) used in the definition of D. From the net number of crossings, 
one can obtain the index of D (|53|). This is the spectral flow method used in 
Refs. []^, ^ to obtain the index of overlap Dirac operator. 



In Fig. 1^, we plot the spectral flow of eight lowest-lying ( near zero ) 
eigenvalues of if^(mo) in the interval 0.9 < mo < 1.3, for one of the gauge 
configurations. Evidently, the flows are not as smooth as one may have ex- 
pected. In this case, the net crossings is -6 ( from negative to positive ), so 
the index of the overlap Dirac operator is -6. In some cases ( see, for example. 
Fig P ), there are some intriguing eigenvalues lying very close to zero, thus for 
a coarse scan in mo, it may not be so easy to determine whether they actually 
cross zero or not. These ambiguities can only be resolved by tracing them 
closely at a finer resolution in mo. 

The distribution of the indices of D for these 100 gauge configurations is 
listed in Table So the index susceptibility is 



4., {in+-n-)^) 



a\ = ^ ^ ' ' = 3.67(50) x 10'^ (54) 
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/ = n+ — n_ 


number of configurations 


4 


2 


3 


11 


2 


9 


1 


14 





24 


-1 


15 


-2 


11 


-3 


6 


-4 


6 


-5 


1 


-6 


1 



Table 1: The distribution of the indices of overlap Dirac operator for 100 
gauge configurations at = 5.80 on the 8^ x 24 lattice. Here (J) = —0.19 and 
= 4.51. 

where A^^ is the total number of sites. Now if we subtract the mean value 
(n+ — n^) = —0.19 from all indices, then the index susceptibility becomes 

a\ = ^ M = 3.64(50) x 10"^ (55) 

Now substitute the index susceptibility ([55|), the lattice spacing (^Oj), /^^ = 
132 MeV, and Nf = 3, into the exact relation (|49|) , we obtain the t]' mass 

rrirj' = (901 ± 64)MeV , (56) 

which agrees with the theoretical estimate 



m2, + m2 - 2m^ = 853 MeV , (57) 

with experimental values of meson massesQ : m^/ = 958 MeV, = 547 MeV, 
and rriK = 495 MeV. 

Next we substitute (|55D and fT^a = 0.0984(3) into (^, and get 

5 = 0.197 ±0.027 , (58) 

which is in good agreement with the value (|ll|) determined from the pion mass, 
as well as with the theoretical estimate S ~ 0.2. 



^Here we distinguish the physical 77' mass ( m,,/ ), from the 77' mass ( m^y ) in queneched 
chiral perturbation theory. 
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4 Concluding remarks 



In this paper, we have examined the quenched chiral logarithm in lattice QCD 
with overlap Dirac quark. The coefficient of quenched chiral logarithm ( 6 = 
0.203 ± 0.014 ) extracted from the pion mass agrees very well with that ( 6 = 
0.197±0.027 ) obtained from the index susceptibility of overlap Dirac operator. 
Further, they are in good agreement with the theoretical estimate 6 ~ 0.2 
in quenched chiral perturbation theory. This provides strong evidences that 
lattice QCD with overlap Dirac quark realizes quenched QCD chiral dynamics, 
as depicted by quenched chiral perturbation theory. 

Our results on pion mass and decay constant rely very much on the via- 
bility of computing quark propagators to a high accuracy. We find that the 
accuracy in the implementation of the inverse square root of is the most 
crucial step in this two- level conjugate gradient paradigm. The Zolotarev op- 
timal rational approximation together with projection of high and low-lying 
eigenmodes enables us to control the error of sign{H^)Y always less than 10^^^, 
at each iteration of the outer CG loop. The details of our implementation are 
described in Section 2. 

Equation (^) is an exact result following from the definition of the 7]' 
mass in (|i5|). It provides the best way to obtain the rj' mass in quenched 
chiral perturbation theory, as well as to determine the coefficient of quenched 
chiral logarithm in terms of the index susceptibility of any topologically-proper 
Ginsparg- Wilson lattice Dirac operator. Equation (|D may even suggest that 5 
is scale invariant for a range of lattice spacings including the continuum limit 
a — ^> 0. This also explains why the S value we obtained at a = 0.147 fm is so 
close to the theoretical estimate in the continuum. 

At this point, it may be interesting to recall an example of lattice Dirac 



operator with exact chiral symmetry |2g|. Unlike the overlap Dirac operator. 



it does not have topological zero modes for any nontrivial gauge backgrounds. 
However, it reproduces correct axial anomaly in the continuum limit at 
least for the trivial gauge sector. Since its index is zero, its index susceptibility 
must be zero, and it follows that the rj' mass is zero, and the quenched chiral 
logarithm is absent in the pion mass as well as other related quantities. Nev- 
ertheless, it does not necessarily imply that this lattice Dirac operator could 
not be realized in nature. It seems that the real testing ground for this lattice 
Dirac operator is lattice QCD with dynamical quarks, in which all quenched 
pathologies are absent. 

In passing, we briefly outline the essential features of our computing system. 
The platform is a home-made Linux PC cluster with 18 nodes, built with off- 
the-shelf components. Each node consists of one Pentium 4 processor ( 1.6 Ghz 
) with one Gbyte of Rambus, one 40 Gbyte hard disk, and a network card. 
The computational intensive parts ( matrix-vector operations ) of our program 



are written in SSE2 codes |28|, which can double the speed of our earlier 
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pure Fortran codes compiled by Intel Fortran Compiler 5.0 with maximum 
optimizations including the SSE2 option. The performance of our system is 
estimated to be around 20 Gflops. 

This work was supported in part by the National Science Council, ROC, 
under the grant number NSC90-2112-M002-021, and also in part by NCTS. 
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Figure 1: The pion mass square (771,^0)^ versus the bare quark mass mqa. The 
sohd hue is the fit of Eq. (|l|) . 
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Figure 2: The pion decay constant /ttO versus the bare quark mass mqa. The 
soHd hne is the hnear fit. 
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Figure 3: 



the fit of Eq. (|T]) divided by m^a. 



Y (fTiqa) versus the bare quark mass m^a. The sohd fine is 
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(m^a)2/(mqa) - B ( mq a ) 
2.2 1 



2.0 - 




0.8 ] , 1 

0.01 0.1 1 

log{ rriq a ) 

Figure 4: The extraction of quenched chiral logarithm by plotting 
(m^a)^/ {rnqO) — B{mqa) versus log(mqa). From the slope of the fitted straight 
line, the coefficient of quenched chiral logarithm can also be determined to be 
5 = 0.2034 ±0.0140. 
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Figure 5: The spectral flow of 8 lowest-lying eigenvalues of for the 50th 
gauge configuration. There are 6 crossings from negative to positive, so the 
index is equal to -6. 
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Figure 6: The spectral flow of 8 lowest-lying eigenvalues of if^ for the 8th 
gauge configuration. There is only one crossing from positive to negative 
around mo = 0.98, so the index is equal to +1. 



23 



